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Abstract 

By using an algebraic diagonalization method, the XYZ Heisenberg antiferromagnetics un- 
der an external magnetic field is studied in the framework of spin-wave theory. The energy 
eigenstates are shown to be squeezed number states and the energy eigenvalues are obtained 
in some cases. Some quantum properties of the energy eigenstates, and the connection of the 
model with the two-mode coupled harmonic oscillators are also discussed. 
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I. INTRODUCTION 



The Heisenberg model [Heisenberg 1928; Anderson 1951], describing spin lattices with short- 
range interactions, plays an indispensable role in the study of the magnetism of solids, and many 
other related physical fields, such as non-local spin systems, quantum dot, nuclear spin, and so 
on [Loss 1988; Burkard 1999; Imamoglu 1999]. According to the sign of interaction intensity 
J, this model can be classified as ferromagnetic (J > 0) type and antiferromagnetic (J < 0) 
type, the latter is much complicated than the former in the eigenstate structure. Based on the 
interaction intensity along different space directions, the Heisenberg model can be classified 
as XXX, XXZ, XYZ, where the former two cases can be regarded as the special case of the 
latter one. In theoretical study, a basic task is to do the digonalization of the Hamiltonian. 
However, due to the complicacy of the many body problem, getting the exact solution is very 
difficult, and till now only some special cases can be exactly solved, such as one dimension 
XXX antiferromagnetic chain using the Bethe-Ansatz method [Bethe 1930]. In practice, the 
spin- wave theory is widely employed as an acceptable approximation, and the "Neel states" 
is regarded to be the approximate ground states of the antiferromagnetic case. For XXZ 
antiferromagnetic Heisenberg model, spin-wave theory has been well studied in an analytical 
way which can be found in standard textbooks [Callaway 1976]. However, for XYZ 
far as we know, there is no result parallel to XXZ case in literature. The reason lies in that, 
although it's not very difficult in principle to deal with the problem through the traditional 
Bogoliubov-valatin transformation [Bogoliubov 1958; Valatin 1958], the solution form is too 
complicated, to be intuitive in physics, and one have to rely on pure numerical analysis, leading 
to the lack of physical transparency. 

As is well known, algebraic method is very important in studying the system with certain 
dynamical group [Kim 1988; Artoni 1991; Pan 2001, 2004]. In this paper, an algebraic diago- 
nalization method will be employed to do the digonalization of XYZ antiferromagnetics under 
an external magnetic field in spin-wave framework. The energy eigenstates which turn out to 
be squeezed number states [Yuen 1976; Kim 1989; Nieto 1997] and the corresponding eigenval- 
ues in some cases will be obtained. The paper is organized as follows. Firstly, the algebraic 
diagonalization method is briefly introduced in section II. In section III, we study the system 
and give the result of eigenstates and eigenvalues in some cases. The relation between algebraic 
method and Bogoliubov method is also discussed. In section IV, we discuss some statistical 
properties of the eigenstates. In section V, we investigate the relation of the system and the 



2 



two-mode coupled harmonic oscillators. Finally, some concluding remarks are given. 



II. AN ALGEBRAIC DIAGONALIZATION METHOD 

The algebraic diagonalization method introduced here is applicable for a class of rather 
general systems whose Hamiltonian possesses dynamical semisimple Lie-algebra structure, i.e. 

H = ^>#* + ]T(A Q £ Q + X* a E^ a ), (1) 

i a 

where {Hi, E a , E^ a } is the standard Cartan-Weyl basis of a semi-simple Lie algebra g satisfying 
the standard commutation relations. 

Based on the Lie algebra g, we introduce a unitary operator 

W(£) = exp[E(Z a E a -&E- a )], (2) 

in which £ a denote parameter functions to be determined. In fact, operator W(£) is a generalized 
displacement operator in the coset space G/H [Zhang 1990]. Utilizing the Baker-Campbell- 
Hausdorff (B-C-H) formula 

e A Be~ A = B + [A,B] + [A, B}} + ■ ■ •, 

and taking advantages of the standard commutation relations of semisimple Lie algebra g, the 
following identities can be obtained 

W(0- 1 G J W(0 = T,^A0G i , (3) 

i 

where Gj stand for the generators Hi and E± a , £ denotes the set of {£, Q }, 7y(0 are parameter 
functions. This is a central step in the overall treatment of the problem. Eq. (3) can directly 
be rewritten into the following form 

W^y'HWiO = J2 e, >^)Hi + E e, \)E a + e, X)E_ a }. (4) 

i a>0 

The non-Cartan generators E± a in the right-hand side are not diagonal and can be eliminated 
through setting 

H a (Z,e,\) = (5) 

for all parameters a which give constraints to the parameter-functions {£q,} and the parameters 
e, A. Then we get 

HW(0 | Ref >= W(OYlvM,e,X)Hi \ Ref > . (6) 

i 
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Here | Ref > denote the reference states required to be the common eigenstates of all Hi. 
Therefore, J2iVi(^ £ ) \)Hi can be substituted by functions acting on the reference states, i.e., 

HW(0 | Ref >= u{£, e, A, Ref)W(0 \ Ref > . (7) 

This is the eigen-equation of Hamiltonian Eq. (1), thus the algebraic diagonalization have 
been accomplished. State W(£) \ Ref > and uj are the energy eigenstates and eigenvalues of 
corresponding system, respectively. 

Now we briefly illuminate the completeness of the solution. For related fundamental in- 
troduction, please see reference [Zhang 1990, p877-879]. In fact, the state | Ref > and the 
operator W(£) here correspond to the state | A, A > (the common eigenstate of Hi) and the dis- 
placement operator Q in that reference, respectively. So, eigenstates W(£) \ Ref > correspond 
to states | A, Q > in that reference which play the role of complete basis of state space. 

III. SQUEEZED NUMBER EIGENSTATES OF XYZ ANTIFERROMAGNETICS 
UNDER AN EXTERNAL MAGNETIC FIELD 

The Hamiltonian of XYZ Heisenberg antiferromagnetics under an external magnetic field 
B = Be z (along z axis) is described by 

H = -J E (V*S?S? + VySfS] + SIS*) + £ B • S, (J < 0, Vx , Vy > 0). (8) 

<i,j> i 

where the notation < i,j > denote the nearest neighbor bonds. We employ the traditional 
two-sublattice treatment, i.e., the spin directions are upwards for sites on sublattice A and 
downwards for sublattice B, and then apply Holstein-Primakoff transformation [Holstein 1949]: 

SI = — s + a) a, s z h = s — b^b, 

Sl = {2s)^a\l-a^a/2sf^ S ] b = (2s)s(l - tfb/2s)h, (9) 

s- = (sty, s b - = (sty, 

where a\ a ($ , b) can be regarded as the creation and annihilation operators of boson on 
sublattice A (sublattice B). The particle numbers a^a, tfb can't exceed 2s. 

In low temperature and low excitation condition, (a/a), (b^b) << s, so the non-linear in- 
teraction in Hamiltonian Eq. (8) can be reasonable ignored [Kittel 1963]. Based on this "big 



4 



s" approximation, transforming the operators into momentum space, and leaving out the bi- 
quadratic terms, we get 

H = 2zsJ[Ns-lC£H k -l)}, (10) 

7~tk = Ok a k + a^k a -k + ^k^k + fr-k^-k + ^k(ak&Lk + a -k&k + a lc^-k + °L k &k) 

+p k (a k b k + a[b[ + a_ k b_ k + a f - k & -k) + A*( & k & k - a lc a k + &L k &-k - a f _ k a -k), (11) 

with 

v k = ^7k, Pk = 7k = -Ee R » = 2^7. (12) 

Here R is a vector connecting an atom with its nearest neighbor, and the sum runs over the z 
nearest neighbors, k is restricted in the Brillouin zone. 2N is total number of the lattices. 
H k can be expressed as the linear combination of six generators of Lie algebra so(3, 2), i.e., 

H k = E* + ^F 3 k + p k (El + £ k ) + v k (Fl + F k ), (13) 

where the generators take the following forms 

E + = a £K + °-k & -k> F + = a k b\ + a- k b^, 

E k = a k b k + a_ k fe-k , = a£b_ k + a\b k , (14) 

E* = l« + n\ + ni + n\ + 2), F 3 k = \{< + ™-k ~K~ n\). 

The other four generators of so(3, 2) are 

C/J = &t k &£, f/ k = M-k, V; k = a+ k a+, V k = a k a. k . (15) 

The non-vanishing commutation relations of this algebra read 

[E + , £_] = -E 3 , [E 3 , E ± \ = ±E±, [F+, F_] = F 3 , [F 3 , F±] = ±F±, 

[£? 3 , £4] = ±£4, [F 3 , U±] = tU±, [E 3 , V±] = ±V ± , [F 3 , V±] = ±V ± , 

[E±, V T ] = TF T , [F ± , U±] = ±E±, [E ± , F ± ] = T V±, [V + , V-] = -(E 3 + F 3 ), 

[F±, V T ] = T E T , [E ± , U t ] = T F ± , [E ± , F t ] = T U ± , [U + , U-] = -(E 3 - F 3 ). 

One can see that {E^,E k ,E^} forms an so(2, 1) ps su(1, 1) subalgebra, while {F k ,F k ,F k } 
forms an so(3) ~ su(2) one. When ^ = r) y , the system will reduce to the su(l, 1) case [Xie 
2002]. 

Hamiltonian Eq. (13) possesses the form of Eq. (1), then one can do the diagonalization of 
Hamiltonian Eq. (13) through the algebraic diagonalization method introduced in section II. 



Note that the reference state | Ref > is required to be the common eigenstate of and _F 3 fc , 
| Ref > should take the form of number state | n k a ,n\ >. So, the eigenstates W(£) | Ref > are 
revealed to be squeezed number state, which had been studied in quantum optics field [Yuen 
1976; Kim 1989; Nieto 1997]. 

In the calculation process of the diagonalization, we found that, if a general full-parameter 
W(^) is set, the infinite progressions induced by the B-C-H formula can't replaced by analytical 
functions. This lead to huge difficult even for number calculation. For the purpose of studying 
the formal analytical solutions and the general property of the squeezed number eigenstate, we 
set 

W(r,9) = exp{r[cos0£+ + sin0(V r + - U+) - /i.e.]}. (16) 

Utilizing the commutation relations of so(3, 2) and the B-C-H formula, after a lengthy calcula- 
tion, we get 

W\r, 6)H k W(r, 9) = cu k £ 3 + oj f F 3 , (17) 

where 

J% = cosh 2r + p k cos 9 sinh 2r, (18) 
oo F = p(l + 2sin 2 6>sinh 2 r) + v k sin 20 sinh 2 r, (19) 

with the constraint equations 

cos 9 sinh 2r + p k (cosh 2 r — cos 29 sinh 2 r) = 0, (20) 
\x sin 29 sinh 2 r + -Uk(cosh 2 r + cos 2 9 sinh 2 r) = 0, (21) 
H sin 9 sinh 2r + p k sin 29 sinh 2 r — v k cos 9 sinh 2r = 0. (22) 

The role of these equations is to eliminate the non-Cartan generators appearing originally on 
the right-side of Eq. (17). Note that there are 5 parameters in Eqs. (20-22), so one of the 
parameters Pk,Vk,/i should be not independent. Acting both sides of Eq. (17) on reference 
number state |n k ,n^), we get 

H k \R k >= K k n k + c k n k + ^ k )|i? k >, (23) 

where 

\R k >=W(r k ,9 k )\n k a ,n*), (24) 



k uj\ + Jp k uj\ - top 

- o ' U b - o • l^J 



Eq. (23) is the eigen-equation. The squeezed number state W(r k , 6> k ) n%) is the eigenstate 
of the system, and lo\ and lo\ are the energies of the two different magnons respectively. It 
is shown that the magnetic field lift the magnon degeneracy. In the study of the number 
calculation, we find usually their exist several solutions, but only one solution give real energy 
eigenvalue. For instance, in a special case the reasonable solution is, p = 0.1, v = —0.0006, 
H = 0.0383, 9 = 1.2506, r = -0.1636, cu a = 0.5418, cu b = 0.5017. 

Now we discuss the relation between the algebraic diagonalization method used here and 
the traditional Bogoliubov-Valatin transformation. To a certain extend, these two methods are 
equivalent. The Bogoliubov-Valatin unitary transformation is known to acts on single creation 
or annihilation operator. In our case, the Lie algebra generators are composed of the quadratic 
form of the creation and annihilation operators. So, the unitary transformation acted on the 
Lie algebra generators corresponds to the product of two unitary transformations acting on 
single creation or annihilation operator respectively. This fact connects the two methods, and 
leads to that the eigenvalues should be the same and the energy eigenstates can be connected 
through considerable transformation. In fact, the unitary operator W(r, 9) (Eq. 16) provides 
a two-mode squeezing transformation, 

a k = W~ x (r, 9)a k W(r, 9) = a k cosh r + (cos 9b\ — sin 9a\) sinhr, (26) 
p k = W-\r,9)b k W(r,9) = b k cosh r + (cos 9a\ + sin 9b\) sinhr. (27) 

Here a k and (5 k are new bosonic operators. In fact, the anti-transformation of Eq. (26,27) can 
just be regarded as a corresponding Bogoliubov transformation. For briefness, we would not 
write them out here. The squeezed number form is (uj^n^+LV^n\+\uj%)W(r^ 9 k )\n^, n%), where 
the n^, are the eigenvalues of the particle number operators of a\.a k , b\b k . The corresponding 
Bogoliubov form is [u^{a k a k + \) + Up(/3 k f3 k + |)]|0, 0), where the a k and f3 k is the bosonic 
operators of the quasi-particles. It's apparent that the two methods work in different operator 
space, one is for " quasi-particle" , the other is for "real particle". 

IV. SOME STATISTICAL PROPERTIES OF EIGENSTATE 

In this section, we discuss some statistical properties of the eigenstate as squeezed number 
state. The average magnetization of sublattice A and B take the forms 

< 4 > = Ns - J2 < b ib >k, (28) 

k 
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< s z A > = -Ns + J2<a j a > k , (29) 

k 

where, 

1 1 

< b ] b > k = -cosh2r(n k + n~ k + n\ + n h k + 2) + -(1 + 2 sin 2 9 sinh 2 r)(n k + - n k - n~ k ) - 1, 

< a f a > k = icosh2r(n k + n~ k + n k + n^ k + 2) + ^(1 + 2 sin 2 sinh 2 r)(n k + ra~ k - n k - n^ k ) - 1. 

When the particle numbers n k and n k vanish, the eigenstate becomes squeezed vacuum state. 
In this state, 

< tfb > k =< at a > k = cosh2r - 1. (30) 

One can see, even in the squeezed vacuum state, i.e., the ground state, the spin direction of 
different sites on a sub-lattice is not necessarily the same; this property of disordering has 
relations with the zero-point motion of the magnon in fact. Moreover, in the squeezed vacuum 
state, the influence of the external magnetic field upon the agree of average spin reversal of 
the two sub-lattice is inviolable the same; while in other squeezed number states, this property 
becomes true only in the condition that n k + n~ k = n k + n^ k . 

Now we study the other quantum effects from the point view of quantum optics. For simplic- 
ity, hereafter only squeezed vacuum state is mentioned. The second-order correlation functions 
of the two excitation modes can be readily obtained as 

<T =rf W = MV)=2 + sin0cothV, (31 ) 
(a k a k ) 



A W = #Wrr = l + 2cos*cothV (32) 



The Mandel Q parameters [Mandel 1979; 1986] read 

gk = gk = (( Aa k fl k) 2 ) _ ! = sinh 2 r + sin Q cosh 2 r (33) 
(a k a k ) 

Eq. (33) illustrates that in squeezed vacuum state, each of the three statistics, sub-Poissonian 
(Q > 0), Poissonian (Q = 0), super-Poissonian (—1 < Q < 0, nonclassical state) are all possible 
to exist in the two magnon modes, relying on the values of the parameters r and 9, which are 
determined by the parameters -u k , p k , /i. 
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If the parameter I k = \J 9a^ / 9ab^ — 1 < 0, then the Cauchy-Schwartz inequality (CSI) 
[Agarwal 1988] is violated, and the correlation between the two k-magnon modes is nonclassical. 
In squeezed vacuum state, 

^ k 1 + (sin# — 2cos#) coth 2 r 

~ l + 2cos#coth 2 r ' ^ ' 

it can be found that larger squeezing parameter r leads to smaller possibility of the violation of 

CIS, i.e. the achievement of nonclassical correlation. This is not surprising since larger squeeze 

parameter r corresponds to larger particle number, tending to approaching classical case. 



V. CONNECTION WITH THE TWO-DIMENSION COUPLED HARMONIC OSCIL- 
LATORS 

For two different physical systems, if their Hamiltonian can be written into the same com- 
bination of the generators of a Lie algebra completely through different realizations of this Lie 
algebra, it's reasonable to believe that the two systems possess the same eigenvalue structure. 
Now we consider the relation of XYZ antiferromagnetic Heisenberg model under an external 
magnetic field and two-dimension coupled harmonic oscillators. Since the so(3, 2) algebra can 
also be constructed in the form 

E + = 44, E- = aia 2 , E 3 = ^(4«i + 4 a 2 + 1), 

F + = ai4, F_ = a[a 2 , F 3 = ^(4«2 - 4°i), (35) 

If 2 1 2 If 2 1 2 

U+ = -a[ , U- = -ai , V+ = -a 2 , V- = -a 2 . 

Eq. (13) can also be realized as 

H! = — — —a\ai H 7p-4°2 + p{did2 + 44) + v (°>i a 2 + 4 a 2)- (36) 

Hereafter we leave out the index k for all the parameters and operators for convenience. So, for 
every k modes, utilizing the realization aj = (mWjXj +ipj)/ \J2mujj, we can map a Hamiltonian 
of two-mode coupled harmonic oscillators corresponding hamiltonian Ji of Eq. (11), 

W ' = S + -7T mx2 j} + \iUxU2mxxXi + A 2 ^^, (37) 



with 



^(A 1 -A 2 )=p, y^(\ 1 + \ 2 )=v, 1-2^ = 2^-1 = /,. (38) 



Hamiltonian 7i and H! should have the same energy spectrum, i.e., in low excitation case, 
every k modes for XYZ antiferromagnetics in external magnetic field can be mapped onto a 
two-dimension coupled harmonic oscillators with x-x and p-p coupling, with the corresponding 
relation of the parameters Eq. (38). One can see that, the system parameters v and p are 
linked to the coefficients of x-x and p-p coupling; while the role of the parameters \i (external 
magnetic field term) is to lift the energy degeneracy of the corresponding harmonic oscillators. 



VI. CONCLUDING REMARKS 



In this paper, based on the introduction of an algebraic diagonalization method, the 
eigenstates of XYZ antiferromagnetics are revealed to be the number squeezed states in the 
framework of linear spin-wave approximation. Different from the traditional Bogoliubov- 
valatin transformation, we study the problem using the algebraic diagonalization method. 
The characteristic of this method lies in that, firstly, one can do calculation in virtue of 
the commutation relations of Lie algebra instead of the Heisenberg algebra; Secondly, the 
energy eigenstates are revealed to be squeezed number states, which had been studied in 
quantum optics field [Yuen 1976; Kim 1989; Nieto 1997], making it convenient to study the 
physical properties of the system. It is shown that, the presence of magnetic field results in 
the non-degeneracy of the two excited magnon; even in the squeezed vacuum state, the spin 
direction on a sub-lattice is not necessarily the same (this disordering in a certain extent is 
attributed to the zero-point motion of the spin oscillator). Moreover, in the squeezed vacuum 
state, the influence of the external magnetic field upon the agree of average spin reversal of 
the two sub-lattice is absolutely the same. Some statistical property of the squeezed vacuum 
state, the second-order correlation functions, Mandel Q parameters, violation of CIS are also 
discussed. By virtue of the algebra method, we illuminate that, XYZ antiferromagnetics under 
an external magnetic field can be mapped onto a two-dimension coupled harmonic oscillators 
with x-x and p-p coupling. Thus, the properties of the two-dimension coupled harmonic 
oscillators with x-x and p-p coupling, can readily be connected to the XYZ antiferromagnetics 
under an external magnetic field through adjusting parameter. Besides, our method can also be 
used to treat the sub-ferromagnetic case, and the case of existing anisotropic crystal magnetic 
field. The application of our results in related physical fields, as well as the realization and 
employment of squeezed number states in other physical systems, leave for further studies. 



10 



ACKNOWLEDGMENT. This work is supported by the National Science Foundation of China 
(11447103), Education Department of Beijing Province and Beihang University. 



Agarwal, G.S., J. Opt. Soc. Am. B 5, 1940 (1988). 

Anderson, P.W., Phy. Rev. 83 1260 (1951). 

Artoni, M. and Birman, J.L., Phys. Rev. B 44 3736 (1991). 

Bethe, H.A. Z., Physik 61 206 (1930). 

Bogoliubov, N.N., Soviet Phys. JETP 34 7 (1958). 

Burkard, G., Loss. D. and Divencenzo, D. P., Phys. Rev. B 59 2070 (1999). 

Callaway, J., Quantum theory of the solid state, Academic, New Yark (1976). 

Holstein, T. and Primakoff, H., Phy. Rev. 58 1098 (1949). 

Heisenberg, W., Phys Z. 49 619 (1928). 

Imamoglu, A. et. al., Phys. Rev. Lett. 83 4204 (1999). 

Kim, M.S., de Oliveira, F.A.M. and Knight, P.L., Opt. Comm. 12 99 (1989). 

Kim, M.S., de Oliveira, F.A.M. and Knight, P.L., Phys. Rev. A 40 2494 (1989). 

Kim, S.K. and Birman, J.L., Phys. Rev. B 38 4291 (1988). 

Kittel, G, Quantum Theory of Solids, Wiley, New Yark (1963). 

Loss, D. and Divincenzo, D.P, Phys. Rev. A 57 120 (1998). 

Mandel, L., Opt. Lett. 4, 205 (1979); Phys. Scr. T12, 34 (1986). 

Nieto, M.M., Phys. Lett. A 229 135 (1997). 

Pan, F. and Draayer, J. P., J. Phys. A., 34 2637 (2001). 

Pan, F., Gueorguiev, V. G. and Draayer, J. P., Phys. Rev. Lett., 92 112503, (2004). 
Valatin, J.G, Nuovo Cimento 7 (1968). 

Xie, B.H, Zhang, H.B and Chen, J.L., Commun. Theor. Phys. 38 292 (2002). 
Yuen, H.P, Phys. Rev. A, 13 2226 (1976). 

Zhang, W.M., Feng, D.H. and Gilmore, R., Rev. Mod. Phys. 62 867 (1990). 



11 



